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Abstract. We consider the Euler-Voigt equations and the Navier-Stokes- 
Voigt equations, which are obtained by an inviscid a-regularization from the 
corresponding equations. The main result we show is the structural stabil- 
ity of the system in term of the variations of both viscosity of rcgularization 
parameters. 



1. Introduction 

One of the most challenging problems in scientific computing is that of produc- 
ing reliable simulations of turbulent flows. Since the work of Kolmogorov in 1941, 
it is well assessed that there are prohibitive limitations due to the smallest scales 
which are persistent in flows at very high Reynolds number. The limitations (due 
to speed but also to memory capacity of the available most powerful computers) 
in performing direct numerical simulations makes the field particularly challeng- 
ing. Incompressible fluid with constant density are described by the Navier-Stokes 
equations 

(la) dtv + {v-\/)v~- lyAv + Vp^ f, 

(lb) V-w = 0, 

supplemented with initial and boundary conditions, where v{t,x) = {v\^V2tV-^ is 
the velocity field, p{t,x) denotes the pressure, f{t,x) = (/i,/2,/3) is the external 
force, and > is the kinematic viscosity. In the sequel we will consider mainly 
the problem in the space periodic setting. 

One main idea is that of studying averaged or filtered equations, see the math- 
ematical overview in [7]. Many models for the numerical simulation of the Large 
Scales (the only ones which are effectively computable and of relevance for the ap- 
plications) has been proposed. Among these models in the recent years there has 
been a lot of activity, from both the pure and applied mathematicians, around the 
so called "a-models." These models are based on a filtering/smoothing obtained 
with the application of the inverse of the Helmholtz operator 

I-a^A, 

and two main questions arise: a) How to describe the (nonlinear) quadratic term 
(I - a2A)-i(u • V) u in terms of (I — a^A) only (interior closure problem in 
LES); b) The role of boundary conditions supplementing the Helmholtz operator 
and the derived model. Roughly speaking, the effect of applying the inverse of the 
Helmholtz operator is that of getting two more derivatives of the solution under 
control. We also recall that Leray's [28] approach to construct weak solutions of 
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the Navier-Stokes equations by smoothing just the convective velocity is based 
on a very similar idea, with regularization made by convolution. There is a big 
variety in the family of a-approximations to the Navier-Stokes equations and we 
recall in alphabetic order, without the claim of being exhaustive, some of the recent 
publications about the subject [ilSIlOlinilHIIlIiaiiailllinilMlISZ]. Further 
details can be found in the introduction of [9l |33l [37] . We observe that essentially 
all the above methods are based on a sort of regularization by a viscous smoothing 
of the equations, which is reflected in better analytical properties. 

A very promising and new approach, recently introduced by Cao, Lunasin, and 
Titi |10| , is that of the inviscid regularization coming from the Layton-Lewandowski 
(or simplified Bardina) model, which is a zeroth order deconvolution method, when 
the viscosity vanishes. Authors in [ini observed that setting the viscosity = in 
that model gives the following (with u = (I — a'^A)~^v) Euler-Voigt model: 



Curios enough, when the viscosity is reintroduced, this turns out to coincide with a 
model for visco-clastic fluids studied starting from the seventies by Oskolkov [Ml 155] 
and known as the Navier-Stokes- Voigt (sometimes written Voight) model: 



The idea of the inviscid regularization is very interesting for two main reasons: 1) 
this is not a viscous regularization, and the energy behavior is respected in a more 
precise way; 2) This regularization docs not introduce new boundary conditions, 
and gives chances to the study of the problem in bounded domain, without the 
aforementioned usual difficulties coming in theory of Large Eddy Simulations, when 
employed in presence of solid boundaries. Related ideas, based on time relaxation, 
has been introduced in |26j , while the same inviscid regularization has been also used 
in the study of water waves [5] , while recent applications to the quasi-geostrophic 
equations are given in [22] . 

A detailed account of many properties of Voigt equations has been given in [171 
[T8l[23l[29] and, beside giving a motivation based on the role in scientific computing 
for the study of the Navier-Stokes and Euler Voigt models, in this paper we treat 
questions more linked with the general theory of partial differential equations, as 
explained for instance in [4]. In particular, the main results, which are obtained by 
using in a simplified setting the techniques introduced by Beirao da Veiga [11 [2l |3] 
concern the well-poscdness of the equations, which is also relevant to understand 
the stability of the equations with respect to small perturbation of the parameters. 
The results we prove are not surprising and they are strictly linked with the similar 
ones recently proved by Linshiz and Titi [29] . Nevertheless our results show how the 
systems are robust and this gives new support for their employment in numerical 
computations. 

Plan of the paper. In Section|2]we introduce the notation and give some remarks 
on the existence of smooth solutions. In Section [3] we prove the sharp limits as the 
regularization parameter a vanishes, while in Section [4] we study both limits as v 
and a vanish. 



(2a) 
(2b) 



dtu - a'^dtAu + {u ■\/)u + \/p ^ f, 

V - u = 0. 



(3a) 
(3b) 



dfU - a'^dtAu + (u • V) m - i^Au + \/p = f, 

V • u 0. 



ON EULER-VOIGT AND NAVIER-STOKES-VOIGT MODELS 



3 



2. Some preliminary results 

Wc introduce now the notation and give some remarks on the existence of solu- 
tions to the Voigt models. 

2.1. Notation. In the sequel we will use the customary Lebesgue and Sobolcv 
spaces VF'''^ and := W^'^, and for simplicity wc do not distinguish between 
scalar and vector valued functions. Since we will mainly work with periodic bound- 
ary conditions the spaces are made of periodic functions. In the Hilbertian case 
p = 2 we can easily characterize the divergence-free spaces by using Fourier Series 
on the 3D torus: Define J7 := [0,27r]^. We denote by (61,62,63) the orthonormal 
basis of M.^, and by x := {xi,X2,X3) E M."^ the standard point in M.^^. Let T be the 
torus defined by T := R'^/Z"^. We use || • || to denote the L^(T) norm and wc impose 

the zero mean condition / (j)dx ^ on velocity, pressure and external force. We 

Jn 

define, for an exponent s > 0, 

Hs := l^w.r ^R^, w e H'iT)^, V • w = 0, Jwdx = 0^, 

where H'iTf [H'{T)]^ and if < s < 1 the condition V • w = must be 
understood in a weak sense. For w G Hs, wc can expand the velocity field with 
Fourier scries 

w{x) = {/}fc6*'"'^, where k is the wave-number, 

and the Fourier coefficients are given by Wfe = pry Jj i(j(x)e~*'^'^c?a;, where |T| denotes 

the Lebesgue measure of T. If k := \k\ = ^/\ki p + |fc2p + |A:3p, then the Hs norm 
is defined by 

where, as above, ||w|1ho The inner products associated to these norms are 

We finally characterize Hg C iJ^^T) as follows: 

Hs := |w = ^ Wfce''=''= : ^ \k\'^'\wk\'^ < 00, k ■ Wk = 0,w^k = }• 

Moreover, we will use the smoothing operator defined for all (5 > and for all 
functions / periodic and divergence free, as follows: For any 7ti > and given 
/ g Hm-i then the function fs defined as 

(4) fs:^ /fc^'""' 

|fc|<l/5 

is such that fs G Hs for all s € M (hence it is infinitely differentiable) and fs ^ f 
in Hm as (5 — > 0. Other properties of this (truncation in wave-number space) 
smoothing operator can be easily obtained from the Fourier series characterization 
and will be recalled in Section [31 

As a final remark on function spaces, we will also use (but only in Section l2.3p 
the characterization of divergence free subspaces of and i/^, with vanishing 
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Dirichlet boundary conditions for which we refer for instance to Constantin and 
Foias [T2] . 

In the sequel (especially to obtain estimates for solutions in Hs, with s non- 
integer) we will also use some elementary commutator type estimates as the fol- 
lowing lemma concerning the operator s G K+ (see e.g. [20l [211 ISIDj where 
A:= (-A)i/2. 

Lemma 2.1. For s > and 1 <r < oo, and for smooth enough u and v 

W{uv)\\Lr < Ci\\u\\LPi\\A'v\\L^, + \\v\\lp2\\A'u\\l^2), 

where 1/r = 1/pi + l/qi = l/p2 + I/92 md C is a suitable positive constant. 

2.2. Existence results. Concerning the Euler-Voigt equations we have the follow- 
ing result of global existence and uniqueness of solutions proved in |231 Thm. 3.1]: 
Therein the result is proved in two different ways, one with the contraction princi- 
ple and the other one by means of the Galerkin method. The results explain the 
hyperbolic nature of the problem and are a starting point also to prove Gevrey 
regularity results. Here we give some remarks on one main technical point and also 
explain some (non strictly essentials) details on the external force. 

Theorem 2.1. (Larios and Titi [23|) Let T > and let Uq G Hm and f G 

C(-r,T; VK™-i'6/5)^ ^^^f^ V • / = 0, /or TO > 1. Then, there exists a unique so- 
lution of the Euler-Voigt equations (Pa|) - (|2bp which belongs to C^[—T,T];H,n)- 
Moreover, 

Wu'^mHr^. <C{a,\\uo\\H^, sup ||/(i)l|m-1.6/5,r), 

-T<t<T 

for all t G [-T, T] . 

Proof. We do not claim any originality in the result, but we just give an alternate 
proof of one main point, emphasizing also the role of the pressure. We use some 
classical techniques employed also in Beirao da Veiga [1] even if here due to the 
space-periodicity the proof is much simpler. We give the proof only for m = 1, 
because this is the most important step. The higher regularity can be obtained in 
the same way by a bootstrapping argument. 

By using the Galerkin approach with approximate solutions (where is 
made with a finite Fourier expansion) it is elementary to show the a-priori estimate 

CGL°°(-T,r;ifi). 

By using standard tools it immediately follows that when m — > -|-c» the approxi- 
mate functions converge to which is the unique solution solutions to ()2a|)- 
(j2bp and which belongs to L°°{—T,T; Hi). One difficulty consists in passing from 
boundedness to continuity with respect to the time variable. This can be obtained 
in an elementary manner as follows: First observe that the pressure p satisfies the 
Poisson equation 

-A39" = V • [(u" • V) u"] . 

Since u" G L°°{-T;T; Hi) we have that V • [(u" • V) u«] G L°°{-T,T;W-'^'^/^), 
where W~^'^^'^ := (W^-^)' and consequently, by using the classical regularity theory 
for the Poisson equation in the periodic setting with zero mean value, we have 

Vp" G L^{-T,T-L^'^). 



ON EULER-VOIGT AND NAVIER-STOKES-VOIGT MODELS 



5 



Then, by comparison it follows that 

and hence by using again the elliptic regularity to the equation (I — q;^A)u" = F 
and the Sobolcv embedding we have 

< e L°°(-T,T;M/2,6/5-) ^ L°°{-T,T;H^). 

By standard results it follows that can be identified with a function continuous 
with values in H^. Coming back (with the improved regularity on u") to the 
estimates on the convective term and on the gradient of the pressure they are both 
now in C{~T,T; L^^^) and the same argument shows that 

< - a^A^ = -Vp" - (it" • V) + / e C(-T, T; L^'^), 

hence that e C{~T,T;H^). This ends the proof, since V • = 0. The high 
order regularity can be proved following the same approach. □ 

To conclude, we also recall the well-known results for the three-dimensional Euler 
equations 

(5a) dtu + {u-V)u + Vp = /, 

(5b) V • u = 0. 

i.e, when a = 0. It is well known that if uq G Hg, and / G L^{Q,T,Hs) with 
s > 5/2, then there exists a unique solution to these equations in C{[0,T]; Hg) n 
C^{[0,T];Hs-i) for some finite time < T < T (see, e.g., the review in ISI|). Let 
us recall that, contrary to the Euler- Voigt equations for which we know existence 
of global smooth unique solutions, it is an outstanding open problem determin- 
ing whether smooth solutions exist globally in time or blow-up in finite time. In 
particular, the best know criterion for the 3D Euler equations is that of Beale- 
Kato-Majda. In [23j Thm. 5.2] it is obtained an interesting criterion linking the 
regularity of the Euler equations, with the dissipation of the Euler- Voigt equations 
as q; — >■ 0. More precisely it is shown that if 

sup limsupa2||Vu"(0||^ > 0, 

te[OT[ Q->o+ 

then the Euler equations starting with the same initial datum uq of {M"}a>o develop 
a singularity in [0,T]. This result is a by product of a result of convergence as 
a —7- 0^ of the Euler- Voigt equations. Consequently, the behavior as a vanishes is 
relevant also in view of applications of this model to computations for the study of 
the possible blow-up for the Euler system. In fact in [531 Thm. 5.2] it is also proved 
that the solutions u" of the three-dimensional Euler- Voigt equations converge to the 
corresponding solutions u of the three-dimensional Euler equations, with respect to 
the norm L°°(0, T; L^(T)), and with initial data = uq G H'lT), for s > 5/2. Our 
main interest is to study the sharp convergence, that in the same space of the initial 
datum. Here, the situation is a little bit different from the usual "Navier-Stokes 
— >■ Euler limit" , since the regularity requested on the initial conditions changes in 
the two systems: The Euler- Voigt requires one more derivative, in order to have 
uniform estimates in terms of a. As a by product of our results we also treat the 
behavior as the positive viscosity v converges to zero and, when introducing the 
viscous problem, we also prove a result on time-periodic solutions. 
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2.3. Time periodic solutions. In this subsection we give some remarks on the 
proof of existence of time-periodic solutions. The results presented here are ob- 
tained with well-established techniques introduced in Prodi , taking the chance 
also of making some observations on the existence for the Navier-Stokes-Voigt equa- 
tions. Since the results proved here hold true also in a bounded domain with 
Dirichlet boundary conditions, we use -just in this section- the notation 

Hi,„ ^ {w.n^M?, w€ H\nf, V • w = 0, w\9n = O} , 

and := {Hi^^)' . Since we are in a case very similar to the 2D Navier-Stokes 
equations (for the Naiver-Stokes-Voigt equations it is easy to prove existence and 
uniqueness of solutions for all times), we can work directly on the solutions, looking 
for a fixed point argument in the infinite dimensional space Hi. For the Navier- 
Stokes equations, the difference between 2D and 3D (regarding time-periodic so- 
lutions) are explained in [30l Ch. 4], where also the Galerkin method with the 
Brouwer fixed point is used to construct approximations to periodic solutions. 

Proposition 2.1. Let C M"^ 6e a bounded, smooth, and open set; let f e 
L^(0, T; i/_i). Then, there exists at least a solution to the Navier-Stokes-Voigt 
equations 

(6a) dtu-a'^dtAu+{u-V)u-iyAu + \7p = f innx[0,T], 

(6b) V-u==0 mnx[0,T], 

(6c) u = ondnx]0,T], 

such that u{0) ~ u{T) G Hi u. 

Remark 2.1. The same result holds also for f2 = T with periodic boundary con- 
ditions. Moreover, the uniqueness of the time-periodic solutions still represent an 
open problem, posing the same difficulties as those well-known for the 2D Navier- 
Stokes equations. 



Proof. Since the proof is very standard, just a sketch of the proof of Proposition l2.1l 
is presented here, as a remark on what can be proved for the Navier-Stokes- Voigt 
equations. We observe that by the same usual methods based on Galerkin approxi- 
mate functions and Aubin-Lions compactness tool we can construct a weak solution 
such that u e L°°(0, T; i/i^^) with 

||m(T)||W||Vu(T)||2+z. r\\Vu{s)fds<\\uo\fW\\Vuo\f+Cr\\f{s)\\j, Js. 

JQ Jo 

Observe that the proof of this result can be obtained by making use of the compar- 
ison argument to prove that 

Ut-a^Aut e L^{0,T;H_i), 

hence that Ut G L^{0, T; Hi^^) by the standard Lax-Milgram lemma set in the space 
Hi^a. In particular, this proves that the solution u e C{\}^,T];Hi,^). 

In order to prove existence of time-periodic solutions we have show that if 
||m(0)|P -I- a^||VM(0)|p < for a large enough i? > 0, then the same bound 
holds true at t = T. Taking the inner product of ((5a|) with u, and by using the 
Poincare inequality we get 

+ a2||V^.|p)+Cl(z.,^],a)(|l^.|p + a2|lV«f) <C2(^,f^)|l/|lL,- 
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Consequently we have that 
e^^^{\HT)f + a^Vu{T)r) < (h(O)f + Vu(0)f ) + c, I ||/W|lL, dt . 



T 

2 







Therefore, to conclude it is sufficient to impose 



l_e-=ir' 

to show that the solution satisfies 

\\u{T)f + a^\\Vu{T)f < 

The proof follows by observing that the ball B{0,R) C -ffi.o- is a convex set in an 
Hilbert space, and it is therefore compact in the weak topology. Hence, by using 
Tychonov theorem we can argue that there exists a fixed point of the map uq — 
u{T), which is then a T-periodic solution to the Navier-Stokes-Voigt equations. □ 

3. Convergence to the solutions of the Euler equations 

In this section we prove the main result of the paper, that is a precise convergence 
result of smooth solutions of the Euler- Voigt equations to smooth solution of the 
Euler equations. 

We start with the following result which is not optimal since one derivative is 
lost in the convergence. The technical reason, which can be easily understood, is 
that the iJm-estimates for the solution of the Euler- Voigt equation starting from a 
datum in H„i are not independent of a. In fact, in this case only the boundedness 
of Q!^||u"||ff^j is known. To have estimates independent of a (and continuity up to 
t = 0) one needs to assume more regularity on the initial datum. On the other 
hand, convergence in i?m-2 is relatively easier to be obtained and the loss of two 
derivatives can be understood from the presence of the term A.uf. We observe 
that this kind of results are obtained in [53]. New technical difficulties arise in our 
setting, since loosing one derivative is in some sense the best result when the initial 
datum is the same for both Euler and Euler- Voigt equations. 

We first prove an intermediate result, since it represents the main technical point. 
Later on we will elaborate on the results which can be obtained when also the initial 
data can change (especially in terms of their regularity). 



Theorem 3.1. Let u be the solution to the Euler equations (|5ap - (j5bp with initial 
condition uq G H^, and let be the solution to (|2ap - ()2bp with initial condition 
Uq — Uq. Let T > be a common time of existence for both u and u°' , with 
u, u"' € C{[0,T]; H3) n C^{[0,T]; H2)- Then, for any sequence {an}, with a„ > 
and such that a„ — > as — > 00, it holds that 

sup |ju""(t) — — > 0, as n—> 00. 

0<t<T 

Proof. The proof uses the smoothing tool introduced in [31 [3] and explained for the 
Euler equations in the periodic setting in [5], see also [32| . Let uq e 7/3 and define 
uq^s as in (j4]), then V • uq^s = and moreover, by direct computation, 

C C 

WuoMIhs < C\\uo\\h3, \\uo^s\\Hi < -J and \\uq,s\\h5 < 
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In addition, for s such that < s < 3, it also holds that ||uo,i5 — uqWh, < CS^^'^. 
Let us be the solution of the Euler equations with initial condition uo_s, which we 
will call "regularized Euler equations." Then, in the interval [0, T] the following 
relations hold true 

C 

(7) \\us{t)\\H., < C and \\u5it)\\H„, < with m > 3. 
We write 

||it" - u\\h2 < - us\\h2 + \\us - u\\h2 =■■ I + II, 
and we estimate \\u°' — '«5||i=o(o.T;H2) ^^'^ ll""? ^ w||L=°(o,T;ff2)- 

Estimate for — us\\l°°{o,T;H2)- We denote by cuf := us — the difference be- 
tween the solution u" of the Euler- Voigt equations ([2a| - (j2bp and us- For simplicity 
we will use the notation lu := ui'g. Thus, we get 

dtoj - o^dtAij + Vips - p") = -a^dtAus - [us ■V)us + {u" ■ V) 

= -a^dtAus - {uj-V)us- (u" • V) w, 

and by taking the i72-inncr product with cj we obtain 

(8) IjtiMk+^'W^'^Wk) <a^\{Adtus,uj)H2\ + \i{uJ-y)us,u;)H2\ 

Next, we estimate the first term from the right-hand side. Let us recall that us is 
a solution to the (regularized) Euler equations, so it follows that 

a^\{Adtus,uj)H2 \ ^ a'^\{dtus,Auj)H2\ 

us, Auj)h2 I + I (^P; Alo)h2 I 
= a2|(V[(u5 • V)] U5,Vcj)jf,|, 

due to periodicity and also to the incompressibility constraint. By using the regu- 
larity of the solution of the Euler equations, we have then 



a2|(V[K • V)u5],Vw)hJ < Ca^{\\usfH^ + \\us\\H2\\us\\HM'^^\\H2 

<Ca\\\us\\l.^+'^^^)\\VLuU2 



<C-\\yu^\\H2, 

where we are supposing for simplicity that < i5 < 1 (since we will main use 
values of 5 close to zero). By using classical estimates on the convective term as 
in Kato [111 Eq. (2.1)-(2.2')], we estimate the other terms from the right-hand side 
of ([8|) as follows: 

\{{u^ ■y)Lo),u:)H2\<\K\\H.AiofH^ and \{{lo ■^)us),lo)h2\<\\us\\h,M?h2- 

Collecting the above estimates, using the bounds for the solution of the (regularized) 
Euler equations, and with Schwarz inequality we get 
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Thus, by using the Gronwall's lemma we infer that 

\K - MUio,T:H.) <{^T + \\uo - uoAl. + «'I|VK - uo,s)\\k)Ce^^ 

(9) 2 



Estimate for \\us — uj|x,oo(-o y.^^)- Here, we take the _ff2-energy estimate for a;^ 
Us — u and we find (since they are both solutions to the Euler equations with 
different initial data) 

^ ■ V)us,u}s)h2\ + \{{u ■ V)ujs,ujs)h2\ 

<Ci\\u\\H, + \\us\\H,)\\^5fH2- 

that is < \\ujo.s\\h2^'"'^ < SC2{T), and finally 

(10) \\us-u\\L^^o^T;H2)<SC2iT). 

As a consequence of the estimates ©-(ITU]) we can conclude that 

where C{T) = max {Cl^'^ (T), C2{T)} . Now, by taking 6 = (5„ such that both (5„ 
and go to zero as n goes to infinity, we obtain the required convergence in the 
i?2-norm. In particular, it follows that — w|lL~(o,T;i/2) = 0{y/a). □ 



u\\L-io.T;H2) < C{S' + + a'\\V{uo,s - iio)||?^j'^'c(r). 



By using exactly the same arguments one can easily prove, more generally, the 
following result in smoother spaces. 



Theorem 3.2. Let u be the solution to the Euler equations (|5ap - (|5bp with initial 
condition Uq G -ffm+i, with m integer such that m > 2, and let u" be a solution 
to the Euler-Voigt equations (j2ap - (|2bp with initial condition Uq = Uq. Let T > 
be a common time of existence for u and u" , with u^^u € C{[0,T]; Hm+i) H 
C^([0, T]; iJm). Then, for any sequence {q;„} with q;„ > and such that — > 
as n 00, it holds that 

sup |ju""(t) — m(<)||h,„ 0, as n ^ CO. 

0<t<T 

We arrive now to the main result of the paper which shows the precise con- 
vergence in terms of the initial data and of the space without unnatural loss of 
regularity. 



Theorem 3.3. Consider the Euler equations (|5ap - (|5bp with initial condition uq G 
H3, and let T > be a finite time of existence for the solution u G C([0, T]; H^) H 
C^([0, T]; iJ2)- Let u"'^ be a solution to the Euler-Voigt equations (|^a|) - (|2bp . with 
initial condition Uq, such that 

i) G Hi, for I5>0 
(II) 

ii) ||uq — uo||h3 — > as (3^0. 
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Then, for any sequence with /?„ > and such that l3n ^ as n ^ oo there 

exists {an} with an > and converging to zero such that 

sup ||u""''^"(t) — u(t)||H3 — > 0, as n —> oo. 

0<t<T 



To prove the Theorem 13.31 we need the foUowing prehminary lemma showing 
that we can construct the solution m""'^" in a time interval independent of n S N, 
and this would be enough in order to get weak convergence results by using the 
classical compactness methods, even if we are interested in strong convergence. 

Lemma 3.1. Under the hypotheses of Theorem \3.3[ it follows that, for any positive 
sequence such that ^„ — ^ 0, as n oo, we can find a positive sequence {a„} 

with an such that ||u""'^"||ioo(Q y./z^-) results bounded uniformly m n G N. 

Proof. Consider the Euler- Voigt equations (|2a |) -(j2b |) with initial data Uq. Taking 
the iJs-inner product with u"'^ , and with the usual inequalities for the convective 
term we obtain 

Thus, we deduce that 



l~Ct{\\u^J%^+a^VuPjj,f^f^ 

with t £ [0,T]. Thanks to the properties i) and ii) in (|TT1) . it follows that, let- 
ting (3^0, then ||wo||i/3 \\uo\\h3- Next, we can choose a ^ such that 
a||VMQ|j//3 remains bounded (for instance choose = OdjuQU^J). ) in this way 
||M"'^||^oo(o,T;ff3) results to be uniformly bounded as well, in a time interval [0,r] 
independent of /3. □ 



We can now give the proof of Theorem [ 
Proof of Theorem \3.3[ To estimate |ju"''^ — u\\ lc^(^o. t-.Hs), we write 
h"'^ - u\\h, < h"'^ - ^^^11^3 + - "11^3 =■! + II, 

< Wu'^'^ - usWhs + \\us - u^Whs + \\u^ - u\\hs =■■ Ii+h+ II, 

where is the solution of the Euler equations with initial data Ug and ug is the 
solution of the Euler equations starting from the regularized initial datum uq^s- 
(Note that, Lemma [3.11 applies also to , and consequently Hu'^Hhj results to be 
uniformly bounded with respect to /3). 

Estimate for — us\\l'>°{o,T;H3)'- By setting uj := us — u"'^ , we get 

dtLU + a^dtAuj + V{ps - p'^'^) = -a^dtAus - (w • V) - {u"'^ ■ V) lu. 

Taking the i/a-inner product with u), and recalling that is solution to the (reg- 
ularized) Euler equations, we arrive at 

+ a^\{V{us-V)us,Vu;)H,\- 
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From ^19] Eq. (2.2)], we get 

|(K^^.V)c.,C.)H3l<h"'^||ff3ll^llH3- 

With direct computations we obtain 

a'^\{W{us ■ W)us,Wuj)h3\ < Ca'^(\\us\\H3\\us\\H4 + II"<5||h2 ll^allHs)!! ^^Ili/a 

< Ca^S-' + S-')\\Vu\\h, = a^Cid)\\Vuj\\H,, 

and also 

|((w • V)u5,w)h3| < WusWhsWojWh^ + \\u5\\hJ\oj\\l'^\\oj\\hs- 

Now, to estimate the second term on the right-hand in Hg', with | < s' < 2, we 
employ usual techniques and Lemma |2. II to get 

(12) ^1(II^IIh.+"'IIv^IIhJ<(II"^!Ih,,,, + II^^"-''IU,,JI1^I1?^^, 

+ a^\iV[ius-V)us],Vuj)H^,\, 

and again 

a^\{\/[{us ■ V)us],Vuj)h^,\< Ca^i\\us\\H^,J\u5\\Hs + \Ws\\h^,J\u5\\h,)\\'7lo\\h^, 

< Ca^{\\us\\H^,^, + S'-'')\\Wuj\\h^, 
<Ca^{l + S^-'') ^a^C{S)\\Vuj\\H^,. 

Then, it follows 

|(||..||1,^, +a2|lV^|ll,J < Cdlc^lll,^, +a2|lVc.|ll,J + «2C2(,S). 

Thus, from the above differential inequality we get an estimate for ^, and 

taking the L°°-norm on [0,T] we find 

Summarizing the previous estimates, we have 

j^iMl, + a^VcWlJ < C{\\us\\h... + WW^-^WnMlk + a'Ci6)\\Vu;\\H, 

+ ^^(||..o||l,^, +a2(||^o||l^^,^^+C2(^)T))^||..||H3. 
To handle the term jlciJollff^, from the right-hand side we write 

(13) Ikollff,, = \\uq.s - u^Wh,, < ||mo,<5 - "o^allif,, + IIuq ,5 - UoIU.'- 

where Uq g is the regularization of the initial datum Uq. Expanding Hmq ^ — ito.a , 
and \\u^ g — UqWh^i in terms of their Fourier coefficients, we have that 

l<|/c|<i 

Therefore, using (jl3p and the above inequalities, we obtain 

(14) \\u^o\\j,,<C6''''\\\4-^o\\l, + \\4\\lJ- 
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Analogously, we get \\ujo\\h^,^^ < CfJ^^^"' - uqHh, + \\u^\\hJ- Then, the differ- 
ential inequality for the iJs-norm becomes 

j^iMk + c^'W^^Wk) <C{\\us\\h. + \\u"^'\\HM\k 

+ a^C{S)\\VLu\\H,+C'{T,a,(3,d)\\Lu\\H,, 

where 

d{T,a,l3,6) := (^S'-^-^\s' + a'){\\u^, - u,\\j,^ + + 
After some manipulations we get 

and by the Gronwall's inequality it follows that 

whereC(r,a,/3,5) -.^ C{\\us{s)\\ Hs + ||u"'^(s)||_H-3 + 1) ds. Now, we have that 

||Vu;o|||^3 < CWViuls - uoMk + G||V(<5 - 4)\\k- 
Then, expanding ||V(?iQ ^5 — uo.is)!!!;^, in terms of its Fourier coefficients, we find 

i<\k\<i 

\k\'\<u~^oA'<YM-M\k^ 

i<ifci<i 

and it follows that 

(16) \ 

<C^\\u^,-uo\\k+Ca'H\\k- 

Hence, by using jlwojll^^ < C'Huo/ — uoIIhj+CII'^o — "olllfj and (fTB)) . the estimate 
becomes 

(17) < [{a^C\6) + C^{T,a, P,6))T + C\\uo.s - uoWk 

Estimate for \\us — U'^\\[^x(^q x:H3)- Writing the iJa-energy estimate for ug — u^ 

and using, for the sake of brevity, uj instead of ujg we have that 
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and by the usual estimates 

(18) j^M\%, < C{\\us\\hs + \\u^\\hs)MI, + C\\us\\hJH\l^Mh,. 

Next, to estimate the second term on the right-hand side of (|18l) we use again an 
i^s'-energy inequahty, with | < s' < 2. Thus, arguing as in the derivation of (fT2)) 
we get 

d_ 

di' 

and using relation (jl4[) . we infer that 
Thus, relation (ITSl) becomes 



where ^(T, /3) := C{T){\\u'^^\\H., + \\u^Q-UQ\\H3f'^- Therefore, applying the Gronwall 
lemma, we deduce that 

\\us - u''|li«=(o.T;_ff,) <(||uo,d- - WoIIh:, + ll"o ^ ""oIIhs 

(19) A 

+ (52-^'C(T,/3)r)e^(^'^'^), 
with C{T,P,5) := C{\\us{s)\\h, + \\uf{s)\\H, + l)ds. 

Estimate for \\u"'^ — u^\\l°='[o,T:H3)'- 
Using dTT]) and dH]) we have that 

\\u"-'^ — M^|lioo(o,T;ff3) ^ C(||u"''' - W<5||ioo(o,T;H3) + \Ws — \\\^ (Q.T-H^)) 

[a^C^S) + C^{T, a, (3, S))T + C\\uo,s - uo\\''h^ 

2 



< c 

(20) 

-4- f :i \ -i- 



+ C{l + -)\\u^,-u4l^+o?U\\l^ 



^2C{T,l3,5)+C(T,a,/3,S) 



Estimate for \\u^ — u||L=«(o.T:ff3): We split // as follows 

(21) 11 <\\u^ -usWHs + Wus-ullHs- 

Considei the first term on the right hand-side of (PT|) . It follows that the difference 
11^5 — u^|!L°°(o,T;_f/3) can be estimated as done in ([T9l). For the second term on the 
right hand-side of (PT|) . we can use the same [35J Eq. (23)]. Hence, we actually get 

\\us - 'u||l=-(0,T;H3) < (11^0,5 - UqWh:, + 5^^"" T)C{T). 

with I < s' < 2. Now, multiplying the right-hand side of the latter inequality by 
gC(T,^,<5)^ and then adding to ([T^ . we get the estimate 

\\u^ - ■"||l°°(0,T;H3) <c[\\uo.S - Uq\\h3{C{T) + C) + \\u'^ - UqWhs 

+ S^-'' {C{T) + CiT, /3))r) e^^^'^'^^ . 



(22) 
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Note that, the terms C{T),C(T, (3), and C{T,l3,d) are bounded in terms of S 
and 13. Then, for any positive sequence {/?«}, such that ^„ — >■ as n — >■ oo, letting 
(5 — > 0, we obtain 

Iju'^" - u||i«=(o,T;H3) asn-^oo. 
Furthermore, by relation ([20]) . letting a„ 0, and taking a sequence (5„ -> 0, such 
that II u^"^^ - uo\\h3, Q;„||u^"||ff4, ff, 17, and go to zero as n goes to infinity, 
we find 

||ya„,/3„ _ ^;3„ ll^^^g ^0 as n — > 00. 

Hence, we have that 

11"""'''" - tillL^Co.Tiffa) ^ asn^oo, 
and the thesis follows. □ 



4. Convergence of the solutions of the Navier-Stokes-Voigt 

equations 

Combining the results of the previous section with similar computations, we 
study also the behavior of solutions in terms of the viscosity. We are still set in the 
space-periodic case and, for simplicity we assume / = 0. Next, we state a conver- 
gence result for solutions of the Navier-Stokes-Voigt equations to the corresponding 
solutions of the Euler equations. 



Theorem 4.1. Consider the Euler equations (|5a[) - (|5bp with initial condition uq & 
H3, and let T > be a finite time of existence for the solution u G C'([0,T]] H3). 
Let u'^''^''^ be a solution to the Navier-Stokes-Voigt equations ()6a|l -([6c |) . with initial 
datum Uq, satisfying the properties i) and ii) in pT|) . Then, for any choice of 
positive sequences {/3n} and {i^n}, both converging to zero as n — > 00, there exists 
a positive sequence {an} converging to zero as n 00, such that 

sup ll'u""'''"'''" — u||/f., — > 0, as n —> 00. 

0<t<T 

Proof. In the following, u°''^ and u"'^-'^ will indicate the solutions of the Euler- 
Voigt and the Navier-Stokes-Voigt equations (with the same initial datum u^) re- 
spectively. To estimate Hm"'^'" — m||l°°(o,T;H3)j we take into account the following 
terms 

(23) 11""'^^"' - "IU-(0.T;//3) ^^^ll""'^'" - <IU~(0,T;//3) 

+ II"? - "llL-(0,T;if3) =■ 1+ II, 

where Ug is the solution of the Euler- Voigt equations with respect to the regularized 
initial datum uo,5- Then, the term // will be split as follows 

// < IK ~ ""•^IU=o(o,T;i/3) + 11""'^ - "llL-(o,T;ff3) Ih + Ih- 

The above splitting is probably not the simplest one, but it is the most convenient 
to employ the results proved in the previous section. 

Estimate for ||u"''^''^ — "a Ili^^co.Tiffg): Setting uj"''^'" := Ug — u"''^''^ (as usual for 
the difference o;^'^'", we drop the symbols a,/3,S, and i'), we get 

dtuj - a^dtAuj - vl\uj + V(p? - p"'"') = -t^Au^ - (w • V) - (m"'^-'' • V) uj. 
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1 d 



Taking the i/a-inner product of the above relation with w, with the same inequahties 
employed in the previous sections we arrive at 

(24) ~iMk+<^'\\^-\\k) + < ^l|V.?||l,3 

In order to estimate ||i^||l=o(t3), we use again the same tool (with the Hs'-enevgy 
inequality, for | < s' < 2.) By Lemma [53] wc get 

^"^/S."!!,, Ml, .,11? 

< Ci\\uj\\l^, + a''\\Vuj\\l^,) + Ci^. 

Then, it follows that 

MUio,T:H^,) < {i\\^o\\l^,+a^Vu;oVH^,) + CiyT)CiT). 
Now, it holds that for < fc < 3 

ll^oillf, = ||uo,5 - u^Wh, <C\\uo,5 - Wo,5llk + ^W'^o.s - ""o Ilk 

<cs^~'HH-uo\\k + \\4\\k)- 

In particular, for k ^ s' and fc = s' + 1 we get 

ho,.-^^o'llH, <C^'^'"'^'(ll4-"o!|l3 + ll"o!lff3)' 

Consequently 

Ml^io,T-M,) < {S'-'-'iS' + Ca')i\\4 - uof„^ + \\uX„^) + Ci.T)C(T). 

We use the above inequality and the bound ||Vw^||_f/3 < C\\u°'\\h3/ 5 < C/6 for the 
solutions to the Euler-Voigt equations (the proof is similar to that of (O). Inserting 
in relation (j24p gives, after some manipulations. 



±i\\c.\\l^+a'\\W^\\lj < C{M\\h. + h-'^'lks + l)Ml^ +C^ + C\T,a,S,,y) 

<C{a,f3J,:.){M%^+a'-\\\/io\\%^)+C^ + C^T,a,S, v). 

Where C{a,l3,5,v) := \\u1\\h., + ||u"''''''||ff3 + 1 and C{T,a,5,v) {S^-"^-'' [5^ + 

Ca^) (||m(^ - -uo 111^3 + 11?^^ II I^J + ^r)^C3(r). Then, using the Gronwall lemma, 
we deduce that 

ll-lli==(o,T;«3) < ((II-0III3 + "'llv-o||l,3) + (c^ + a^(r,a,<5,.))r)e^(^.-''.*», 

with C{T,a,p,5,v) := C i^{\\u1{s)\\H, + ||u"''3^'^(s)||h3 + l)ds. 
To conclude, we need a further estimate of q;^|| VwoHl/., = a^||V(uo,(5 — ^o)lllf3- 
Consider the regularized initial data ^. It follows that 



(26) 



. „,/3 l|2 , /^„2|i„,/3 „,/3|i 

Hi 



a2||V(uo,5 - 4)!Ih3 < Ca^\W, - <,||1,^^ + Ca1|<, - 4f 

<c-^\\4-u4i,+co?\\<rH,. 
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and noting that \\uq^s - Uq\\%^ < C\\uo^s - uqWh^ + C\\uq - uoW'j^^, we obtain 

^^^^ lhlli~(0,T;H3) <C(ll"0,5 - UoWh.^ + C(l + -^)\\uo - Uo\\h^ 



52 

Estimate for — u°'''^\\l^(^q^t-H3)' Taking the i/a-energy estimate for ' "'^ 
u" — u"'^ (also in this case, we drop a, (3 and 6), we arrive at 
d 

We estimate the term ||w||^oc(t3) in the usual way and we find 

Then, applying the Gronwall lemma and using the bound (j25p . we get 

making computations similar to those performed in (jl8p - (|19p . we finally obtain the 
estimate 
d 

dt' 



where d{T,a,(3,S) {S^-^'' {6^ + a^){\\4 - ^oWji, + cHt). Hence, we 

get 

Now, arguing as in (p6)) . we can conclude that 



'2 

^2^^ IK - ^"'^llioo(o,T;H3) <c{\\uo,s - uoWl, + C{1 + ^)\\u^o - uoWl, 
where C(T, a, /3, (5) ~ C j'^ {\\u'I{s)\\h, + \\u"^P{s)\\h, + l)ds. 

Estimate for Wu""'^ — u||ioo(g y.^^-): Note that, up to a sub-sequence a„ — >■ as 
n — > 00, the needed estimate on ||it""'^'> — m||loc(o,T:_H3) is provided by Theorem l3.3l 
As a consequence of the above bound and the estimates (P?]) and (pS)) . for any 
positive sequence {vn\, with — >■ as n — >■ 00, letting /?„ -> 0, we can choose 
a pair of sequences q;„,(5„ (look at the proof of Theorem 13. 3p . such that 
ll'«o,5„ - uo||h3> "rJl^o'llffi' 17' -^rrr, and go to zero as n goes to infinity. 
Hence, we get 

||^a„,;3„,i/„ _ — > as n — > 00, 
and the thesis follows. □ 
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Remark 4.1. The result holds also in spaces of more regular functions Hm, m > 3, 
hy using essentially the same techniques. 
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